Generalized stacking fault energetics and dislocation properties: compact vs. spread 

unit dislocation structures in TiAl and CuAu. 



00 
OS 



Oleg N. Mryasov 1 , Yu. N. Gornostyrev 1,2 and A.J. Freeman 1 
1 Department of Physics and Astronomy, Northwestern University Evanston, IL 60208-3112. 
2 Institute of Metal Physics, Ekaterinburg, Russia. 

We present a general scheme for analyzing the structure and mobility of dislocations based on 
solutions of the Peierls-Nabarro model with a two component displacement field and restoring forces 
determined from the ab-initio generalized stacking fault energetics (ie., the so-called 7-surface). 
The approach is used to investigate dislocations in Llo TiAl and CuAu; predicted differences in 
the unit dislocation properties are explicitly related with features of the 7-surface geometry. A 
unified description of compact, spread and split dislocation cores is provided with an important 
characteristic "dissociation path" revealed by this highly tractable scheme. 
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I. INTRODUCTION 

In the development of modern materials [With ad- 
vanced high-temperature structural propertiesoo, there 
is a growing awareness of the importance of a fundamen- 
tal understanding of the role of specific lattice defects. 
One such defect, dislocations, is always present in most 
real materials and influences a variety of properties, in- 
cluding mechanical, electric, optical and magnetic. High- 
temperature intermetallics is one of the most vivid exam- 
ples where particular dislocations control deformation, 
fracture and other properties that are most important in 
terms of their technological applications. 

Dislocations are a unique example of the extended de- 
fect for which a single length scale description (atomistic 
or continuum) is not quite appropriate - a feature that is 
one of the main sources of difficulties in their theoretical 
description. There are hundreds of atoms in a dislocation 
core (requiring an atomistic treatment) and thousands 
of atoms are involved in the formation of the elastic dis- 
lacements fields (where a continuum description is more 
appropriate) so that neither entirely atomistic nor contin- 
uum descriptions are suitable to the problem but rather 
need to be combined. 

Classical atomistic simulations, however, proved to be 
very useful and provided important insights into their 
fundamental properties in metalsH. To make use of atom- 
istic simulations for modeling dislocations in particular 
material reliable model of interatomic interaction has to 
be developed. This, however, poses major principal prob- 
lem for the use of classical atomistic simulations. Unfor- 
tunately, in most of interesting cases (including inter- 
metallics like Llo TiAl) interatomic interactions are very 
complex and so their satisfactory description can not be 
guaranteed withing existing central force models. This 
significantly limits abilities of classical atomistic simula- 
tions for predictive analysis dislocations properties par- 
ticularly in intermetallic alloys with their complex char- 
acter of interatomic interactions!]. To date, the few quan- 



tum molecular dynamic simulations of dislocations have 
been restricted to s,p systems (Si, SiC), employing peri- 
odic boundary conditions with a few hundred atoms in 
the supercelfi. This approach remains unexceptably ex- 
pensive and has principal difficulties to be applied for 
intermetallics containing 3d-elements. 

The complexjty of the problem has stimulated re- 
newed interestETQ in the simple but tractable Peierls- 
Nabarro (PN) dislocation model, which incorporates 
atomistic effects, in an approximate way into a contin- 
uum frameworka. In this model, the effects of lattice dis- 
creteness are confined to a slip plane with a local relation 
for non-linear atomistic restoring forces which can be de- 
termined from an effective rigid shift interplanar poten- 
tial known as a generalized stacking fault (GSF) energyu. 
The PN model proved to be very useful and has strongly 
influenced conceptually the theory of dislocations, frac- 
ture and related mechanical properties! Naturaly ap- 
peared in PN model, the energy of the relative displace- 
ments of two halves of the crystal (ie., the so-called gen- 
eralized stacking fault (GSF) energy, or 7-surface) was 
introduced and studied as an important general charac- 
teristic of the mechanical response in solidaapa .— ■ 

Despite its well-known approxmrntionsEJallStirl^i the 
PN model is very attractive since it offers a tractable de- 
scription with parameters (7-surface) accessible by highly 
accurate quantum mechanical calculationsEm. As was 
recently demonstrated, the applicability of the Peierls 
model framework can be extended in a semi-discrete ver- 
sion even for such an extreme case as narrow dislocation 
cores in Sili3. In this work, the classic Peierls model dis- 
location energy functional is minimized using a discrete 
scheme which allows avoiding assumptions of the dislo- 
cation core planarity and the Nabarro approximation for 
calculating Peierls stresses. Unfortunately, the proposed 
semi-discrete scheme is purely numerical; however, in the 
theory of such complex lattice defects as dislocations it 
is highly desirable to keep the high tractability of the 
theoretical approach as much as possible. 
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To date, the PN model analysis with first-principles 
parametrization for GSF energetics was restricted, how- 
ever, to considering dislocations with a single component 
displacement field (one dimensional (1D1 PN model ) in 
SB and in B2 intermetallic compounds!! The abilities 
of the PN model to deal with dislocation cores in inter- 
metallics has not yet been studied systematically, except 
for our earlier overall encouraging experience with the 
1D-PN model for dislocations in B2 NiAl and FeAlQ. 
However, the two component displacement field (two- 
dimensional (2D) PN model) with the possibility of dis- 
location dissociation is quite typical situation in general 
and, in particular, is necessary to be considered in the 
case of unit dislocations for Llo TiAl and CuAu, ma- 
terials with similar fcc-like structure but very different 
mechanical behaviour. 

In contrast with the ID PN model, exact solutions of 
the 2D-PN model are unknown even for the simplest si- 
nusoidal restoring force law. Recently, Schoeck presented 
a direct variational method for analyzing, on the energy 
functional level, the 2D-PN-podel using ansatz solutions 
with geometric parameters!!]. This method has been used 
to calculate dislocation core structures corresponding to 
model GSF surfacesEH. Here we present an alternative 
scheme for solving 2D-PN model equations within a wide 
class of analytic functions which makes this scheme uni- 
fied (i.e., free of parameters and assumptions on disloca- 
tion structure), suitable for analyzing realistic/complex 
7-surfaces and physically transparent. We apply this 
scheme to analyze unit dislocation structure and mobil- 
ity in Ll TiAl and CuAu alloys with GSF energetics 
determined from ab-initio calculations. 

In this work, we demonstrate how advantages of the 
PN model's high tractability can be fully revealed due to 
the use of a general and physically transparent scheme for 
analyzing solutions in the generalized case of a two com- 
ponent displacement field (2D) and with a first-principles 
determination of the 7-surface energetics. We focus on an 
analysis of the unit dislocation structure in Llo TiAl and 
CuAu using the proposed scheme since this type of dis- 
location is important for understanding these materials 
with very different mechanical behavioura. We demon- 
strate that this proposed approach explains qualitative 
differences in dislocation structure and mobility in a very 
natural way and establishes an explicit relation between 
7-surface geometry and dislocation properties within the 
2D PN model. This allows one to relate GSF energet- 
ics determined from ab-initio electronic structure calcula- 
tions with dislocation structure using a unified and highly 
transparent procedure. 



II. APPROACH AND METHOD 

The structure of dislocations within the 2D PN model 
is described from a balance between elastic and atomistic 
restoring forces which is expressed in a system of integro- 



differential equationalall^ 
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Here /3 is the index for components of the dislacement 
u, withp£t=l for screw and /3=2 for edge, and K a p = 
[iDp6 a piB, Dp is a parameter equal to 1 for (3 = 1 and 
1/(1-^) for /3=2. There are at least two difficulties to 
applying the PN model for realistic situations: (i) the 
restoring force law is unknown in general and (ii) an 
analysis of Eq. (|l|) seems to be difficult for an arbitrary 
restoring force law. We follow the GSF conceptfl and. 
define the restoring force in the local approximation£L3 
as Fp(u) = —dE(u)/dup, where GSF energies E(u) are 
calculated for experimental lattice structure parameters 
employing the local density full-potential linear jauffin- 
tin orbital (FLMTO) band structure technique^!] with 
the Ceperly- Alder form for the exchange-correlation po- 
tential. 

In fact, special points on the 7-surface known as pla- 
nar stacking fault energies (APB, ISF, CISF etc.) have 
been calculated by various ab-initio methods for a num- 
ber of metals, alloys and compounds (see references int3). 
Thus, calculations of the entire 7-surface do not pose 
principal difficulties but may be rather challenging from 
a numerical stand point if one would like to achieve, as 
much as possible, equally high accuracy for all 7-surface 
sampling points. Most of the numerical aspects one may 
face in ab-initio calculations p£ the GSF energetics have 
been summarized by Paxtoncj. Probably, the first ex- 
ample of the entire 7-surface determined using ab-initio 
techniques is a pseudopotential plane wapj (PPW) cal- 
culation for Si by Kaxiras and DuesberyE-3. The PPW 
approach, however, is somewhat less natural to use if 
rather localized d and f electons are present in a sys- 
tem. Thus, in this work we make use of the all-electron 
FLMTO technique which allows one to treat equally well 
both localized and delocalized electronic states-with rea- 
sonable accuracy and high numerical efficiency!!!]. 

One should emphasize here that only in the local ap- 
proximation is so simple a relation provided with GSF 
energies and hence restoring forces are so easily accessible 
by ab-initio techniques. As demonstrated by Miller and 
Phillipsli3 using MD simulations, the local approximation 
merely affects dislocation structure but may be impor- 
tant for calculating Peierls stresses. Any corrections be- 
yond the local approximation involve a much more com- 
plex proceduret^l and to the best of our knowledge have 
not yet been realized using ab-initio calculations. 

An analysis of the solutions of Eq. (pi) is rather complex 
from a mathematical point of view. An elegant idea pro- 
posed by Lejeck for finding solutions of the ID PM model 



equations in a wide class of analytic functional!! seems 
to be general and suitable for cases of complex/realistic 
restoring force laws. Details of such an analysis for the 
ID PN model with generalized restoring forces deter- 
mined for NiAl from ab-initio calculations are discussed 
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in our earlier workS. Our results for this realistic example 
demonstrate that the scheme for solving the PN model 
type equations based on Lejeck's idea can be very conve- 
nient and provides the accuracy and numerical stability 
desired. In the 2D PN model case, however, additional 
complications related with the coupling of equations (cf. 
Eq.(l])) have to be overcome. In this work, we present a 
scheme for solving the 2D PN model equations with an 
arbitrary restoring force law which is based on Lejcek's 
ideaEil as generalized for the 2D case. 

Solutions of Eq. ([!]) are found by constructing an an- 
alytic complex function Gp(x) = pp(x) + i ■ ■jjyFpfa) 
where pp(x) — dup(x)/dx is the dislocation density, and 
Fp(u(x)) the restoring force and are given as expansions 
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Here N is the number of poles with maximal order at 
points in the complex plane, z% — + «Cf (k = 1, N). 
This representation is not only mathematically justified 
(the structure of the PN model equations is analogous 
to the Hilbert transform) but is also physically transpar- 
ent. Indeed, the number of poles in the expansion for 
pp(x) determines the number of cores (partial or frac- 
tional dislocations) with a width £^ and positions of the 
centers, x^. The solutions, up(x), can be found by in- 
tegrating pp{x) with appropriate boundary conditions. 
Then, for an fee lattice and a screw orientation of the 
unit dislocation in the {111} plane, the functions up{x) 
and Fp(u{x)) for p&=2 and N=2 can be written in the 
convenient parametric form 
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Here 6>f = 2arccot((x — x%)/(P), C, 13 = apu^, where luq 
is the width of the dislocation core which appears in the 



solution of the original PN model for a sinusoidal restor- 
ing force law, b is the total Burgers vector and b e is the 
edge component of the partial dislocation Burgers vec- 
tor. Thus, the complete set of dislocation structure pa- 
rameters, a (} and x^ (including possible splitting^ can be 
determined from the general restoring force lawE3, Fp(u). 

However, there is a complication related with an es- 
sential feature of the 2D PN model. Now, the 2D PN 
model equations (Eq.(|])) are coupled onlyta through the 
restoring force which is a function of the unknown a pri- 
ori "dissociation path" determined in the two compo- 
nent displacement space and given by the dependence 
u 2 = f{ u i) °f the edge (u 2 ) on the screw (m) compo- 
nent. Thus, function f(u\) has to be found along with 
solutions of Eq.(|l|). Hence, we propose an iterative proce- 
dure which starts with calculations of the restoring forces, 
Fp(ui,U2), and the corresponding dislocation structure 
parameters for a linear "dissociation path"cJ, u 2 = K-u\. 
Then, as convergence is achieved, both the actual " disso- 
ciation path" and the solutions u(x) can be determined. 



III. RESULTS AND DISCUSSION 

The results of such an analysis for a unit dislocation 
structure in Llo TiAl and CuAu are presented in Fig. 
|l|. As can be seen, a unit dislocation in TiAl is compact 
- in contrast with CuAu. The unit dislocation struc- 
ture predicted here for TiAl (comnact core) is consistent 
with availablfi-|experimental datao and atomistic simu- 
lation resultscJ (for an embedded atom potential that 
gives a value of the CSF energy that is similar to our ab- 
initio calculations) but our analysis is performed within a 
highly tractable scheme which allows one to relate quali- 
tative differences in dislocation properties to the ab-initio 
7-surface geometry. 

Dislocation dissociation and ^-surface geometry - Fo- 
cus now on the details of the 7-surface geometry which 
determine the difference in the structure of TiAl and 
CuAu unit dislocations. These 7-surfaces, calculated 
for the {111} plane with the FLMTO method using 
six layer supexcaUs and homogeneous, periodic bound- 
ary conditionscj'Lj without relaxationEH, are presented in 
Fig. |^ as contour plots along with actual "dissociation" 
paths. Note that the dislocation axis is taken along (110) 
so that displacements along (112) correspond to an edge 
and those along (110) have a pure screw component. We 
label several special points on the 7-surface section along 
(112) for the screw component ui=0.5: (i) the complex 
stacking fault (CSF), labeled C, (ii) the inflection point, 
labeled I, and (iii) the point C'=(0.5,f(0.5)), where f(0.5) 
is the calculated maximal edge component. There is a 
noticeable difference in the relative location of I and C 
for TiAl and CuAu which results in striking differences 
in the properties of the unit dislocations. 

If the components of the restoring forces, Eq. (g), 
along the u 2 = f(ui) path change sign then there is 
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more than one pole in the expansion, Eq. (g), and so 
solutions of Eq. (|l|) correspond to the partial/fractional 
(split/spread ) dislocations. In particular, this condi- 
tion is fulfilled if the corresponding 7-surface section 
has a minimum at (ui=0.5, U2=f(0.5)), where the maxi- 
mal edge component £(0-5) was found to be substantially 
smaller than the trial2J already after the first step of the 
iterative procedure. This "reduction" is a result of the 
overlap of the opposite sign partials displacement fields. 
As a consequence of such a "reduction", C may be lo- 
cated either before the I point (CuAu) or beyond (TiAl); 
this, in turn, determines whether there is a minimum on 
the 7-surface section and, correspondingly, whether the 
dislocation is split/spread or compact. 




-1 



1 



FIG. 
tion: a 





x, a 

1. Density of displacements p(x) in the unit disloca- 



the screw component in TiAl and for b) the screw and 
c) edge components in CuAu as a function of the distance (x) 
from the dislocation center (given in lattice constant units). 

The iterative procedure allows us to also find the actual 
"dissociation path", curve ui = /(wi), which is shown in 
Fig. H as a thick solid line. In the case of a spread 
unit dislocation (CuAu), the actual .path is quite differ- 
ent from the usually assumed pathcil. However, for a 
compact dislocation (TiAl) the difference of the actual 
"dissociation path" (solid line) from the assumed path 
(dashed line) is really dramatic (see Fig. |J). The possi- 
bility of such a transition from the split/spread to com- 
pact solution (with a corresponding transformation of the 
"dissociation path") has also been recently indicated by 
Schoeck in a 2D PN model analysis with model type 7- 
surfaces using the direct variational approactO. 

The equilibrium separation distance between partials 
of dissociated dislocations is traditionally determined 
from the balance between elastic (within linear elasticity 
theory) and stacking fault surface teasipa forces resulting 
in the widely used simple expressionllSlEl. The PN equa- 



tions (Eq. (|l])) do not contain, in an explicit form, the 
stacking fault energy (7/) which is related in a very sim- 
ple manner with the equilibrium separation-distance (d) 
between partials of dissociated dislocationsEB. However, 
it is possible to demonstrate that this simple relation can 
be recovered within our more general PN model approach 
as the limiting case of large separation distances. 




-0.3 -0.2 -0.1 -0.0 

<112] 

FIG. 2. Contour plots of the 7-surface (in J/m 2 ) for Llo 
a) TiAl and b) CuAu. The position of the CSF is denoted 
by C, the point corresponding to the calculated amplitude 
of the edge component by C and the inflection point of the 
corresponding section by I. The actual "dissociation path" 
(or dependence of the longitudional on lateral to the Burgers 
vector components of the the displacement) is presented by 
the thick solid lines and is compared with the usually assumed 
path (straight dashed line). 

Indeed, the energy 7^ of the fault with vector /=(0.5, 
U2) is uniquely determined within the local approxima- 
tion through the two component restoring force vector 
F(u) defined for an arbitrary path u\ = f(u2) as 



4 



7/ = J F{u)du = J (Fipi + F2P2) dx 



(8) 



(0) 



where 7/ is uniquely defined relative to the reference 
point E(0)=0. The integral in Eq. (||) can be evalu- 
ated analytically using solutions of the PN model in the 
form of Eqs. (Q) and (||), and so gives a one to one corre- 
spondence between 7/ and parameters of the dislocation 
structure, j(uf) = j(Ci,Xi). In the limit of large sepa- 
ration distances (d/xi <C 1), the following two simplifi- 
cations are well justified: (i) the relation between edge 
and screw components is close to linear ( the dependence 
ui = f(v,2) is a straight line) and (ii) the maximal edge 
component is given by the position of the CSF. Then, a 
power series expansion of the general relation (Eq. (Q)) 
in Q/xi (up to second order terms ) has the relatively 
simple form 
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where the separation distance in the limit considered is 
taken as d — 2x\ = 2x2 with a ~ a\ — 012, b s = 6/2 is 
the screw component and b e the edge component of the 
Burgers vector and 0*1,2, (1.2 are parameters in Eqs. @ 
and (Q) . Then it is easy to see that the zeroth order terms 
of the 7/ = 7(C, d) power expansion in £/d (for large 
separation distances 2x\ = 2x 2 = d and ( <^ d) repro- 
duce exactly the well-known relation between d and 7/li3. 
This simple estimate, however, becomes obviously inac- 
curate for dislocations with small separation distances, 
as in the case of CuAu and TiAl considered here. This 
result is non-trivial since it allows one to relate rigoc-i 
ously the PN model and often used simple elasticitvt2l 
descriptions. We also find that the higher order (in (/d) 
contributions in Eq. ^| have negative sign, thus the sim- 
ple relation 7 ~ X/d corresponds to an upper estimate 
which becomes increasingly inaccurate for dislocations 
with small separation distances, d. We should empha- 
size, however, that in accordance with our 2D PN analy- 
sis the CSF energy alone is not a sufficient characteristic 
to describe important details of the dislocation structure. 
Hence, additional parameters of the entire 7-surface (see 
Eq. |9|) have to be taken into account for a correct analysis 
- which is especially important in the case of small dislo- 
cation spreading (as found in our analysis for 7-surfaces 
calculated for TiAl and CuAu). 

The Peierls stress. - Consider now how such a charac- 
teristic of the dislocation mobility as the Peierls stress, 
op, can be determined within the proposed PN model 
analysis scheme. According to the approximations in- 
troduced by Nabarrco (which often work reasonably 
welfj), up is defined as 



d<S> P {l) 
dl 
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Here h is the lattice .period in a direction perpendicular to 
the dislocation line,E3, and <5 is the relative displacement 
of the atomic arrays in neighbouring {111} planes of the 
fee lattice. Hence, the determination of a p is reduced to 
the integration of the product of analytic vector functions 
(Eqs. |) and (1)). 

We used these expressions to calculate Peierls stresses 
for the unit dislocations in TiAl and CuAu and found for 
the screw orientation that the ratio o v j \i (here p is the 
(9)shear modulus) for CuAu (0.0005) is significantly smaller 
than for TiAl (0.024). This difference is attributed to the 
spreading of the unit dislocation in CuAu, since <J p /pi for 
the unspread core solution is about the same as for TiAl. 
This result is quite remarkable since it clearly demon- 
strates the importance of the entire 7-surface geometrical 
parameters as an integrated characteristic of the inter- 
atomic interactions of solids (including chemical bonding 
contributions) in the context of the analysis of dislocation 
properties. 

We emphasize here the importance of such character- 
istics that are naturally revealed in our 2D PN model 
analysis as a " dissociation path" for the analysis of the 
fracture behaviour. As is well-known in the scope of 
Rice-Thomson (R-T) criteria, crack propagation is de- 
scribed by the competition of the crack opening (char- 
acterized by the surface energy 7 S ) and dislocation emis- 
sion (characterized by unstable stacking fault energy 7 US ) 
processes. An analysis of this competing process arrives 
at .thp .criterion of brittle crack propagation, ~f s /jus > 
2.9EirEil, which despite its simplicity captureSj-the gen- 
eral tendencies in materials fracture behaviourEil. Dislo- 
cation emission from the crack tip within the R-T cri- 
teria is characterized by the energy barrier, j us . This 
energy barrier is usually-jestimated assuming a straight 
line "dissociation" pathej. As we demonstrated above, 
the actual "dissociation" path may deviate significantly 
from an assumed straight line. In cases of spread dislo- 
cations (CuAu), corresponding changes in 7„ s estimates 
are likely to be small but may be quite significant in cases 
of compact dislocation cores (TiAl). In fact, if features of 
the actual " dissocation" path are taken into account to 
determine 7 US for TiAl, the R-T criterion {^ s /lus = 1-7) 
predicts brittle crack propagation - opposite to the re- 
sult which one obtains assuming a straight " dissocation" 
path {^s/lus = 4.5). 



through the so-called "misfit energy" per unit length of 
dislocation, §p(l), which can be written in the formQ 
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IV. CONCLUSIONS 

We proposed a general and physically transparent 
scheme for analyzing the 2D PN model with a gen- 
eral restoring force law. We provided realistic examples 
(TiAl, CuAu) and demonstrated that dislocation prop- 
erties can be predicted and explained on the basis of ab- 
initio calculations of the 7-surfaces. In particular, the ac- 
tual "dissociation path" was determined and proved to be 
an essential feature of the 2D PN model and an important 
characteristic for understanding dislocation dissociation 
processes and related mechanical properties (tendencies 
in fracture behaviour). We showed that differences in 
dislocation properties (compact vs. spread) for realis- 
tic (complex) situations as in TiAl and CuAu can be 
explained from the ab-initio 7-surface geometry. This 
demonstrates the advantages of the proposed approach 
which establishes, within the 2D PN model, explicit re- 
lations between 7-surface geometry parameters (which 
can be accurately determined using ab-initio techniques) 
and dislocation structure and mobility (processes that 
approach mesoscopic scales). 
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